Tuning the properties of complex transparent conducting oxides: role of crystal 
symmetry, chemical composition and carrier generation 
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The electronic properties of single- and multi-cation transparent conducting oxides (TCOs) are in- 
vestigated using first-principles density functional approach. A detailed comparison of the electronic 
band structure of stoichiometric and oxygen deficient In2C>3, a- and /3-Ga2 03, rock salt and wurtzite 
ZnO, and layered InGaZnCU reveals the role of the following factors which govern the transport and 
optical properties of these TCO materials: (i) the crystal symmetry of the oxides, including both 
the oxygen coordination and the long-range structural anisotropy; (ii) the electronic configuration 
of the cation(s), specifically, the type of orbital(s) - s, p or d - which form the conduction band; and 
(iii) the strength of the hybridization between the cation's states and the p-states of the neighboring 
oxygen atoms. The results not only explain the experimentally observed trends in the electrical 
conductivity in the single-cation TCO, but also demonstrate that multicomponent oxides may offer 
a way to overcome the electron localization bottleneck which limits the charge transport in wide- 
bandgap main-group metal oxides. Further, the advantages of aliovalent substitutional doping - an 
alternative route to generate carriers in a TCO host - are outlined based on the electronic band 
structure calculations of Sn, Ga, Ti and Zr-doped InGaZnO.4. We show that the transition metal 
dopants offer a possibility to improve conductivity without compromising the optical transmittance. 



I. Introduction 

Multicomponent transparent conducting oxides 
(TCOs) - complex oxides which contain several types 
of main group metal ions - have been developed [lHy] 
to broaden the range of transparent conducting and 
semiconducting materials required for a variety of 
applications, including photovoltaic cells, flat panel 
displays, and flexible and invisible electronics. Binary 
and ternary compounds and solid solutions with elec- 
trical, optical and mechanical properties controlled via 
chemical composition, have been a subject of numerous 
investigations, e.g., see 043 for reviews. Since the 
1990s, the multi-cation TCOs which include metal ions 
beyond the traditionally employed Sir, Cd, In and Zn 
have emerged, for example, Mgln204 [7|, Galn03 [8| and 
the so-called 2-3-3 or 3-3-4 systems where the numbers 
correspond to divalent, trivalent and tetravalent cations 
Q . Among them, layered compounds of the homologous 
series (In,Ga)203(Zn O) ra , w here n=integer, attract an 
increasing attention [10h22i |. originally, due to a possi- 
bility to enhance conductivity via spatial separation of 
the carrier donors (dopants) located in insulating layers 
and the conducting layers which transfer the carriers 
effectively, i.e., without charge scattering. However, the 
electrical conductivity and carrier mobility observed 
to-date in these complex materials, cr=100-400 S/cm 
and ^=10-20 cm 2 /V-s, respectively 0, El EE USEEI 
are considerably lower than those achieved in the 
single-cation TCOs, such as 1x1203, Sn02 or ZnO, with 
cr=10 3 -10 4 S/cm and ^=50-100 cm 2 /V-s. 

The low conductivity in (In,Ga)203(ZnO)„ which de- 
creases continuously as the number of ZnO layers, n, in- 
creases, along with the observed strong anisotropy of the 
transport properties in these complex oxides [13J, LL5j was 



assumed to correlate with the density of the octahedrally 
coordinated In atoms that changes with n. Since late 
1970's, the octahedral oxygen coordination of cations was 
believed to be essential for a good transparent conduc- 
tor 0, IE [E [2E H3 ■ Accordingly, it has been suggested 
that in (In^Ga^O^ZnO^ the charge is transferred only 
within the octahedrally coordinated In layers while the 
Ga-Zn layers where both cations have tetrahedral oxygen 
coordination, were believed to be non-conducting 10i.l28|. 
However, the above assumption appears to be in conflict 
with: 

(i) experimental observations that the mobility and 
conductivity are independent of the compositional vari- 
ations in the amorphous InGa03(ZnO)„ and the conclu- 
sion that this material is a Zn 4s conductor 251 ] . 

(ii) electronic band structure investigations which 
showed that all chemically and structurally distinct lay- 
ers in InGaZn04 (n=l) and in similar compounds are 
expected to participate in the charge transport owing to 
the hybrid nature of the bottom of the conduction band 
formed from the s-states of all cations in the cell and the 
p-states of their oxygen neighbors [ifjj ]. 

The above discrepancy calls for further microscopic 
analysis of the layered oxides. In particular, the underly- 
ing mechanisms associated with carrier generation in the 
complex TCO materials should be addressed. The pur- 
pose of this work is not only to determine the properties 
of oxygen deficient or substitutionally doped InGaZn04 
but also to understand how the presence of several cations 
of different valence, ionic size, electronic configuration 
and oxygen coordination affects the carrier generation 
and the overall electronic and optical properties of com- 
plex TCOs. 

For the materials of the (In,Ga)203(ZnO)„ family, it 
has been accepted that free carriers appear due to the for- 
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mation of oxygen vacancies. Also, it was assumed that 
each vacancy donates two extra electrons, i.e., there is no 
charge compensation via natural defects or defect com- 
plexes. However, it was shown theoretically [3(j that 
uncompensated oxygen vacancy in InGaZnC>4 results in 
a deep level within the band gap and, thus, free carriers 
cannot be efficiently generated. Therefore, oxygen va- 
cancy cannot explain the observed conducting behavior 
in this material. 

Substitutional doping, an alternative route to gener- 
ate carriers, remains unexplored in InGaZn04 and the 
compounds of the homologous series. One of the main 
reasons for this is the structural and compositional com- 
plexity of multicomponent oxides. Indeed, the presence 
of several cations with different valence (e.g., In 3+ and 
Ga 3+ vs Zn 2+ ) and oxygen coordination (octahedral for 
In vs tetrahedral for Ga and Zn) may facilitate a non- 
uniform distribution of impurities including their cluster- 
ing or secondary phase formation. Furthermore, targeted 
doping required for extra carrier generation via aliovalent 
substitution, becomes difficult as the number of the struc- 
turally and/or chemically distinct cations in the complex 
material increases. The aliovalent substitutional dopants 
known to be efficient carrier donors in the constituent 
single-cation TCO, e.g., Al 3+ in ZnO, may not be suit- 
able for the complex oxides since they may prefer to sub- 
stitute the cations of the same valence (In 3+ or Ga 3+ ). 
In this case, extra carriers are not generated. Formation 
of charge compensating defects and defect complexes is 
also more likely in complex oxides. Thus, the range of 
dopants which are efficient for a given multicomponent 
material may be quite narrow. 

In this work, first-principles density functional ap- 
proach is employed to investigate possible carrier gen- 
eration mechanisms in the layered InGaZn04 (n=l) and 
its constituent single-cation oxides, In203, Ga203 and 
ZnO, and to understand how each of the mechanisms 
affects the resulting electronic and optical properties. 
First, we perform electronic band structure calculations 
of oxygen deficient materials to understand the origin of 
the observed anisotropic transport properties in the lay- 
ered multicomponent oxide. For this, we analyze the role 
of crystal symmetry and metal-oxygen hybridization by 
comparing the properties of high- and low-symmetry ox- 
ide phases. In addition, our comparative studies of sev- 
eral oxides allow us to understand how chemical compo- 
sition (In vs Zn vs Ga oxides) affects the resulting optical 
and electronic properties. 

Next, in search of efficient substitutional dopants, 
we study the structural and electronic properties of 
InGaZn04 doped with either conventional dopants such 
as the post-transition metal ions (Sn, Ga) or with the 
transition metals (Ti, Zr). We find that the optical 
and transport properties are both sensitive to the ionic 
radius of the dopant and its electronic configuration. 
In particular, we show that the presence of the local- 



ized d-states is beneficial for balancing good conductiv- 
ity with low optical absorption. Based on the results 
obtained, the (dis) advantages of both carrier generation 
mechanisms - oxygen reduction and aliovalent substi- 
tutional doping - are discussed. We conclude that the 
conductivity currently achieved in the oxygen deficient 
(In,Ga)203(ZnO)„ can be significantly improved with 
proper doping of the host material. 

The paper is organized as follows. First, theoretical ap- 
proaches and structural peculiarities of the compounds 
under investigation are outlined in Section II. Studies 
of oxygen deficient In 2 03, a- and /?-Ga203, cubic and 
hexagonal ZnO are presented in Section III. Here we dis- 
cuss the formation energies of the oxygen vacancy in var- 
ious charge states and the changes in the host electronic 
band structure caused by the presence of the oxygen de- 
fect. Significantly, the range of materials investigated al- 
lows us to address the long-standing question about the 
role of the crystal symmetry, both local (oxygen coor- 
dination) and long-range (structural anisotropy), in the 
charge transport. In addition, wc compare the electronic 
band structure features of the three oxides and draw 
general conclusions about the role of the cation-anion 
orbital hybridization in the resulting transport proper- 
ties. In Section IV, we perform similar investigations for 
the layered undoped stoichiometric and oxygen deficient 
InGaZn0 4 . 

In Section V, aliovalent substitutional doping is consid- 
ered. For each of the dopants studied, i.e., Sn, Ga, Ti or 
Zr, several distinct site locations in the layered InGaZn04 
are compared based on the total energy calculations. The 
results, along with the considerations of the dopant and 
the substituted cation valences, allow us to judge the effi- 
ciency of each dopant, i.e., its capability to donate extra 
electrons. Further in this section, the role of the dopant 
electronic configuration in both the transport and opti- 
cal properties is discussed and practical ways to optimize 
the former without compromising the latter are outlined. 
Finally, calculations for InGaZn04 with interstitial oxy- 
gen atoms located within either the In-0 or Ga-Zn-0 
layers are performed to investigate a possibility of charge 
compensation in substitutionally doped materials. 

We discuss the implications of the obtained results in 
Section VI. 

II. Approach 

First-principles full-potential linearized augmented 
plane wave method }3lL l32j | with the local density ap- 
proximation (LDA) is employed for the electronic band 
structure investigations. Cutoffs for the basis functions, 
16.0 Ry, and potential representation, 81.0 Ry, and ex- 
pansion in terms of spherical harmonics with £ < 8 inside 
the muffin-tin spheres were used. The muffin-tin radii are 
2.3 for In; 1.85 a.u. for Ga; 1.9 for Zn; and 1.6 a.u. for O 
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atoms. Summations over the Brillouin zone were carried 
out using at least 23 special k points in the irreducible 
wedge. For every structure investigated, the internal po- 
sitions of atoms were optimized via the total energy and 
atomic forces minimization. 

To study the electronic properties of undoped stoi- 
chiometric oxides, we also employed the self-consistent 
screened-exchange LDA (sX-LDA) method 33j-|37( for 
more accurate description of the band gap values and the 
valence/conduction band edges. The sX-LDA is known 
to provide a better energy functional beyond LDA or 
generalized gradient approximation (GGA) by modeling 
the exchange-correlation hole within the nonlocal density 
scheme. 

InGaZn04 has rhombohedral Rim crystal structure 
of YbFe2C>4 type [H, |3{J. Indium atoms reside in the 
3(a) positions (Yb), while both Ga 3+ and Zn 2+ atoms 
are in the 6(c) positions (Fe) and are distributed ran- 
domly _J40( - as confirmed by our total energy calcula- 
tions (29(. Because of the different ionic radii, the Ga 
and Zn atoms have different z component of the internal 
site position (0,0,z): we find z=0.214 for Ga and z=0.218 
for Zn. The optimized cation-anion distances are 2.16- 
2.22 A for In-O; 1.88 A and 1.96-2.16 A for the Ga-0 in 
the a6-plane and along the [0001] direction, respectively; 
and 1.97 A and 1.99-2.36 A for the planar and apical 
Zn-0 distances, respectively. These values correlate with 
the ionic radii of the cations (0.94 A for six-coordinated 
In, 0.61 A for four-coordinated Ga and 0.74 A for four- 
coordinated Zn) and correspond to the cation-anion dis- 
tances in the single-cation oxides, i.e., In203 (2.13-2.25 
A), Ga 2 3 (1.83-2.07 A) and ZnO (1.97-2.14 A). 

To model isolated point defects, the following super- 
cells were constructed: a 49-atom supercell with the lat- 
tice vectors (302), (112) and (021), given in the units 
of the rhombohedral primitive cell vectors, and a 77- 
atom supercell with the lattice vectors (302), (131) and 
(212) for InGaZn.04; a 40-atom and 80-atom supercells 
for bixbyite hr^Os; a 120-atom supercell for monoclinic 
/3-Ga2 03; a 80-atom supercell for rhombohedral (corun- 
dum) a-Ga203; a 84-atom supercell with the lattice vec- 
tors (221), (131) and (231), given in the units of the 
hexagonal primitive cell vectors for wurtzite ZnO; and 
64 and 128- atom supercells for cubic ZnO. 

To compare the properties of oxides, we choose the su- 
percells which result in a similar defect concentration, 
namely, 0.8-1. Ox 10 21 cm~ 3 . This corresponds to the dis- 
tance between the oxygen defects of 8-10 A. We stress 
that such high concentrations may not be experimentally 
feasible in all materials under consideration. However, 
our goal is (i) to compare the properties in oxides un- 
der the same (perhaps, hypothetical) conditions, and (ii) 
to demonstrate that even when the defects arc relatively 
close to each other in our model and so their wavefunc- 
tions are expected to overlap significantly and, hence, 
form a broad band, we obtain a localized state (as, for 



example, in /3-Ga 2 03, see below). 

The oxygen vacancy with three charge states, i.e., neu- 
tral Vq and ionized Vq or Vq + , was modeled using a 
corresponding background charge. The defect formation 
energy is calculated as: 



AH(Ep, [l) — Edefect — E, 



host 



;±(m° + a m ) + 5 £ f , 
(i) 

where Edefect and Ehost are the total energies for the 
oxygen deficient oxide and the stoichiometric oxide in 
the same-size supercell, respectively; A/i is the chem- 
ical potential for atom added to (-) or removed from 
(+) the lattice; q is the defect charge state and Ep 
is the Fermi energy. The chemical potentials A/i are 
taken with respect to the LDA energy fjP of the elemen- 
tary metals or the O2 molecule. We apply the thermo- 
dynamic stability condition for the oxides. Our calcu- 
lated formation enthalpies, AHj, are in a good agree- 
ment with the experimental values given in parenthesis: 
AH f {In 2 3 )=-8.2 cV (-9.6 eV); AH f {ZnO)=-3A 
eV (-3.6 eV); AH f (Ga 2 3 )=-9.7 eV (-11.2 eV); and 
AH f(InGaZnOi)—~ 11.97 eV (experimental data not 
available). The formation energies of substitutionally 
doped In203 and InGaZn04 were calculated with respect 
to the bulk orthorhombic Ga, tetragonal In and Sn, and 
hexagonal Ti, Zr and Zn metals. 

III. Electronic properties of I112O3, G&20>, and ZnO 

III.l Undoped stoichiometric TCO hosts 

Oxides of post-transition metals with (n — l)d 10 ns 2 
electronic configuration, such as 1x1203, Ga203, ZnO and 
CdO share similar chemical, structural and electronic 
properties. They have densely-packed structures with 
four- or six-coordinate metal ions, Table HI Strong inter- 
actions between the oxygen 2p and metal ns orbitals give 
rise to electronic band structures which are qualitatively 
similar among all of these oxides, Figure [TJ the bonding 
and nonbonding O 2p states form the valence band while 
the conduction band arises from the antibonding Ms-Op 
interactions where M stands for metal ion. The empty 
p-states of the cation (dashed line) form the following 
band at a higher energy. The partial density of states 
plots, Figure [TJ reveal that the oxygen 2p (thin line) and 
metal ns (thick line) states make similar contributions to 
the conduction band. This provides a three-dimensional 
Ms-Op network for charge transport once extra carriers 



fill the band [41 1 



Ms-Op interactions result in a gap between the va- 
lence and the conduction bands. In ZnO, the gap is 
direct whereas in I112O3 and Ga203 the valence band 
maximum is at the H point ([111]) and M point ([111]), 
respectively, giving rise to indirect band gap of 2.88 eV 
and 4.86 eV, respectively. (For results on CdO, see Refs. 
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FIG. 1. Electronic band structure and density of states (DOS) of undoped stoichiometric (a) bixbyite In203, (b) wurtzite ZnO, 
(c) monoclinic Ga203, and (d) layered InGaZnC>4 as calculated within the self-consistent sX-LDA. The thick, dashed and thin 
lines in the DOS plots represent the metal s, metal p and oxygen p states, respectively. 



In 2 3 ZnO Ga 2 3 InGaZn0 4 



Lattice Bixbyite Wurtzite Rocksalt Monoclinic Corundum Rhombohedral 

Coordination of cation 6 4 6 6,4 6 6(In),4(Zn,Ga) 

Coordination of anion 4 4 6 4,3 4 4 

Fundamental gap, eV 2.88 (1.16) 3.41 (0.81) 2.49 (1.03) 4.86 (2.32) 5.16 (2.81) 3.29 (1.19) 

Optical (direct) band gap, eV 3.38 (1.52) 3.41 (0.81) 3.86 (2.28) 4.91 (2.32) 5.40 (3.06) 3.40 (1.28) 

Electron effective mass, in m e 0.28 (0.18) 0.35 (0.19) 0.28 (0.22) 0.34 (0.24) 0.32 (0.23) 0.34 (0.21) 

Effective mass anisotropy, 5 1.00 (1.00) 1.01 (1.25) 1.00 (1.00) 1.10 (1.23) 1.00 (1.00) 1.00 (0.97) 

TABLE I. Structural and electronic properties of TCO hosts - the undoped stoichiometric oxides. The values are determined 
within the self-consistent screened-exchange local-density approximation and within LDA (given in parentheses). Anisotropy 
of the electron effective mass is defined as 8 = (m lwo] + m [010, )/2m [001, . 



42M45| .) For TCO applications, direct optical band gaps 
are of primary importance. The calculated values are 
given in Table fl] As one can see, our sX-LDA results are 
in excellent agreement with the reported experimental 
and theoretical values, namely, 2.9 eV for the fundamen- 
tal ba nd gap and 3.5-3.7 eV for th e op tical band gap for 
ln 2 3 E^pl, 3.1-3.6 eV for ZnO (Ej and 4.5-4.9 eV for 
Ga^O^ [5lL l52j|. while LDA significantly underestimates 
the band gap values, as expected. 

The Ms- Op overlap also determines the energy disper- 
sion of the conduction band in these materials. Within 
the framework of k-p theory [53|, the electron effective 
mass can be found within the second-order perturbation: 



,( c ) 



E 

l=£c 



(uv\pi\uw)\' 



(2) 



where p is the momentum operator, \u^) is the Bloch 
wave function of the I's band at the T point (wave vec- 



tor k=0) and is its energy. Band label c represents 
the conduction band, while the sum runs over all other 
bands. In the oxides under consideration here, the elec- 
tron effective mass is less than the mass of the electron, 
m e . As it follows from Eq. [2j it is determined primar- 
ily by the valence band contributions (2?W < E^), i.e., 
by the oxygen 2p-states. Thus, the network of alternat- 
ing metal and oxygen atoms ensures the small electron 
effective mass in the TCO hosts 41 1. 

Moreover, the electron effective mass varies insignifi- 
cantly with oxygen coordination which may be different 
for different phases of the same-cation oxide [29]. For 
example, for ZnO in rock salt (octahedral coordination) 
or wurtzite (tetrahedral coordination) structures, and for 
In2 03 in Ia3 (bixbyite), R3c (corundum) or 12 1 3 struc- 
tures, the effective masses vary by less than or about 15% 
(see also [HI). In addition, the effective mass remains 
nearly isotropic in all oxide phases considered including 
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those with well-defined long-range structural anisotropy. 
This is in accord with the investigations of complex ox- 
ides of main group metals 2{J 54 1 where isotropic electron 
effective mass is found for oxides with irregular atomic ar- 
rangements or large structural voids. These observations 
explain the success of (i) amorphous TCOs whose elec- 
trical properties remain similar to those in the crystalline 
state, and (ii) TCO-based flexible transparent conduct- 
ing coatings. 

As shown in the next section where we consider the 
conversion of the undoped stoichiometric TCO hosts 
from insulators into conductors, the Ms-Op origin of the 
conduction band plays critical role in determining the 
resulting transport properties. 



III. 2 Oxygen deficient InzOz, ZnO and GC12O3 

Oxides of post-transition metals exhibit a relatively 
low free energy of formation [55| which favors large oxy- 
gen deficiencies even under equilibrium growth condi- 
tions. This gives rise to the free-carrier densities as large 
as 10 17 -10 19 cm" 3 in ln 2 3 and ZnO [HcliEl]. 

Removal of an oxygen atom from the metal oxide leaves 
two extra electrons in the crystal. Whether one or both 
of these electrons become mobile (free) carriers or remain 
localized at the vacancy site is determined by the forma- 
tion energy of the oxygen defect which can have various 



charge states (i.e., V G , Vq or V Q + ). We find that in all 
oxides investigated in this work, the neutral oxygen de- 
fect Vq corresponds to the ground state of the oxygen va- 
cancy, Fig. [5J This is in accord with the results for 1^03 
and ZnO by Lany and Zunger 58[. From the electronic 
band structure calculations for oxygen deficient In2 03, 
Ga203 and ZnO, we find that the uncompensated oxygen 
vacancy, Vq, gives rise to a fully occupied (by the two 
vacancy- induced electrons), thus, non-conductive state, 
as the Fermi level falls into a gap between the defect 
state and the rest of the conduction band. 

The observed conductivity in L12O3 and ZnO was at- 
tributed to the photoexcited oxygen vacancy, i-e^, Vq~, 
which provides a metastable conductive state [58|. In- 
deed, when the oxygen vacancy is singly ionized (Vq), the 
lowest single conduction band associated with the pres- 
ence of the oxygen defect becomes half occupied, leading 
to a conducting behavior Within the local density 
approximation (LDA), we estimate the optical excitation 



V° 
v o 



V+ + e" to be 1.64 eV, 1.48 eV and 2.02 eV for 



L12O3, ZnO and Ga203, respectively. These values are 
underestimated because LDA underestimates the band 
gap by 2-3 eV, cf., Table H 



Thus, only if the oxygen vacancy is excited to V, 



conducting behavior may occur. The conductive state 
associated with Vq is found to be located closer to the 
conduction band minimum (CBM) as compared to Vq, 
Table [ill Again, the values in Table HI] for the location 



of the defect states with respect to band edges may be 
underestimated, however, we expect that LDA captures 
the trend among the oxides 60] and below we use the 
LDA results to compare with those obtained for oxygen 
deficient a-Ga203, rock salt ZnO and InGaZn04. 

Despite similar defect concentration and, hence, sim- 
ilar distance between the oxygen defects (8-10 A) m 
L12O3, ZnO and /3-Ga203, the obtained widths of the 
lowest conduction band in the conductive oxides (i.e., 
those with Vq defect) differ significantly, Table HU The 
lower energy dispersion of the conduction band in f3- 
Ga203 suggests that the extra electrons tend to local- 
ize - in marked contrast to the extended state in 1^03. 
Accordingly, the calculated Fermi electron group veloc- 
ities in the conductive oxides differ by almost an order 
of magnitude, Table |TTJ We stress that these differences 
cannot be accounted by the variation in the electron ef- 
fective masses of stoichiometric materials (cf., Table [J). 
The obtained trend in the electron group velocities is in 
accord with experimental results: indium oxide is known 
to be the best TCO in terms of the minimum resistivity 
achieved, whereas /3-Ga203 has not been considered as 
a viable TCO and, moreover, is believed to be the con- 
stituent detrimental for the charge transport of a multi- 
component oxide. To understand the microscopic origin 
of the different electronic properties of the oxygen defi- 
cient oxides, further comparative analysis will follow (see 
Sections III. 3-5). 

At the end of this section, we would like to point out 
the differences in the electronic band structures of oxy- 
gen deficient and substitutionally doped oxides. Figure[3] 
illustrates how the electronic band structure of the In203 
host is affected by the presence of oxygen vacancy (Vq) 
or substitutional Sn atom (Sn/„). For the same concen- 
tration and the same charge state of the two defects, we 
find that oxygen vacancy leads to a non-rigid-band shift 
of the Fermi level, and a large second gap appears due to 
the splitting of the lowest single conduction band from 
the rest of the conduction band, Fig. [3] Conversely, 
there is no second gap in the case of Sn substitution. We 
have found earlier [6l| that the second gap previously re- 
ported for Sn-doped In2 03 [62| vanishes upon structural 
relaxation around Sn ions. In marked contrast to the 
substitutional doping, the second gap increases by about 
23 % when the structure of the oxygen deficient 1^03 
is allowed to fully relax. Indeed, the absence of a sec- 
ond gap in case of the substitutional doping indicates a 
good hybridization of the dopant states (Sn s-states) with 
the states of the neighboring atoms (O p-states), while 
an oxygen vacancy implies a substantial redistribution 
of the electron density to accommodate the defect and 
to minimize its impact. The obtained differences in the 
structure of the conduction band bottom of the doped 
and oxygen deficient oxide will, clearly, give rise to dif- 
ferent transport and optical properties. 
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Fermi level E F - E v (eV) Fermi level E F - E v (eV) 

FIG. 2. Calculated formation energies of the oxygen vacancy as a function of the Fermi level. The symbols indicate the 
transition energies between different charge states of the oxygen defect. Solid lines stand for In203 (triangles), wurtzite ZnO 
(diamonds), monoclinic /3-Ga203 (squares) and InGaZn04 (circles). Dashed lines are for rocksalt ZnO (open diamonds) and 
rhombohedral a-Ga203 (open squares). The oxygen defect concentration in all systems is 0.8-l.lxl0 21 cm -3 . The results are 
obtained within LDA. 




FIG. 3. Electronic band structure of (a) undoped stoichiometric In203, (b) 6.25 at.% Sn-doped L12O3 which corresponds to Sn 
concentration of 1.93x 10 21 cm -3 , and (c) oxygen deficient L12O3 with the oxygen vacancy Vj concentration of 1.96x 10 21 cm -3 
as obtained within self-consistent sX-LDA. 



III. 3 Electron localization in oxygen deficient TCO 

Our goal is to understand the microscopic origin of 
the different electronic properties (cf. Table |TTJ) of the 
three oxides in which the cations have the same elec- 
tronic configuration, i.e., s°p°. In particular, we would 
like to understand the origin of the strong electron lo- 
calization that manifests itself in the narrow conduction 
band in oxygen deficient ^-G^O^. For this, we analyze 
the charge density distribution within the lowest con- 
duction band and also compare the contributions to the 
conduction band from the cations and anions located at 



different distances from the oxygen defect. 

For I112O3 we find that the In atoms nearest to the 
vacancy site give 2-3 times larger contributions than the 
rest of the In atoms in the cell. As a result, there is a 
notable build-up of the charge density near the vacancy, 
namely, 31 % of the total charge at the bottom of the 
conduction band reside on the nearest neighbor In and 
the next nearest neighbor O atoms. 

In oxygen deficient ZnO, there is an order of magni- 
tude difference between the contributions from the spe- 
cific atoms (either cations or anions), and 39 % of the 
total charge density belong to the nearest Zn and next 
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TABLE II. Electronic properties of oxygen deficient In203, wurtzite and rock-salt ZnO, a-Ga203 and /3-Ga203, and InGaZnCU 
as calculated within LDA. The conduction band width reflects the degree of electron localization. The smaller the s-orbital 
contribution on the cations nearest to the oxygen vacancy, the more charge stays trapped near the defect. 





L12O3 


h-ZnO /S-GaaOs 


c-ZnO 


a-Ga203 


InGaZn04 


Oxygen defect concentration, xl0 21 cm -3 


0.97 


1.00 


0.80 


0.80 


1.30 


1.11 
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nearest O neighbors of the oxygen vacancy. 

In Ga2C<3 we find that 61 % of the total charge is at 
the defect's nearest neighbor cations and the next near- 
est neighbor oxygen atoms, with the rest of the cations 
or anions giving up to two orders of magnitude smaller 
contributions. 

Thus, the electron density distribution is more uni- 
form within the cell in indium and zinc oxides, whereas 
a strong electron localization is observed in the oxygen 
deficient Ga2C<3. For further analysis, we consider the 
role of the crystal symmetry (Section III. 4) and the type 
of orbitals (s, p or d) which form the conduction band 
(Section III.5). 

III. 4 Role of the crystal symmetry 

The trend in the amount of the build-up charge near 
the vacancy site and the resulting energy dispersion of 
the conduction band and associated electron group ve- 
locities (Tabic [TTJ in the single-cation oxides correlate 
with their crystal symmetry, i.e., cubic In203 > hexag- 
onal ZnO > monoclinic Ga2 03. A higher symmetry is 
expected to provide a larger orbital overlap and, hence, 
a better charge transport. We wish to shed some light 
onto this correlation. 

111.4. A hong-range symmetry: crystal anisotropy 

First, let us consider the role of the long-range crystal 
lattice anisotropy which is exists in /3-Ga2C>3 (monoclinic 
phase, space group C2/m) due to the presence of struc- 
turally distinct (non-equivalent) Ga and O atoms. As 



one can see from Table U the effective mass anisotropy 
is the largest for the undoped stoichiometric /3-Ga 2 03, 
with the masses found to be equal to m[ 100 ]=m[° 10 l=0.35 
m e and m[° 01 ]=0.32 m e within sx-LDA. However, this ef- 
fective mass difference in the oxide host cannot explain 
the obtained strong anisotropy of the electronic prop- 
erties that appears in the oxygen deficient /3-Ga203 as 
discussed below. 

Due to the presence of the non-equivalent Ga and O 
atoms, the distribution of the oxygen vacancies is not 
uniform. We find that Vq prefers to be located in the 
oxygen sites of the 0(1) type, while the total energy of 
the configuration with the defect in the 0(2) or 0(3) sites 
is higher by 163 meV or 267 meV, respectively. The oxy- 
gen atoms of the 0(1) type form zig-zag chains along the 
y direction with the atomic distance within the chain of 
3.04 A, while the shortest distances between the chains 
are 3.91 A and 4.76 A in the x and z directions, re- 
spectively. Therefore, the interaction between the oxy- 
gen vacancies distributed among the 0(1) sites should be 
stronger within the chains than between the chains. This 
leads to different Fermi electron group velocities along 
the main crystallographic directions: we obtain 0.8 xlO 5 
m/s, 1.9xl0 5 m/s and 1.3xl0 5 m/s along the x, y and 
z directions, respectively. The conductivity which is pro- 
portional to the square of the electron velocity is expected 
to be highly anisotropic - with the factor of 6 between 
the x and y directions. 

Thus, we argue that crystal lattice anisotropy affects 
the transport properties of a TCO primarily due to a 
non-uniform distribution of carrier donors whereas the 
electronic properties of the TCO host remain nearly 
isotropic. For each compound, there may be a particu- 
lar carrier generation mechanism found which maintains 
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FIG. 4. Octahedral coordination of oxygen atom by cations 
provides the largest overlap between the s-orbitals of the six 
metal ions and p x , p y and p z orbitals of the oxygen atom. 

isotropic nature of the oxide host. 

III.4-B Local symmetry: oxygen coordination 

As discussed in Section III.l, the conduction band in 
the oxides under investigation is formed from the s or- 
bitals of the cation and p orbitals of the oxygen atoms. 
These orbitals will participate in the charge transport 
once the materials are degenerately doped. The orbital 
overlap between the neighboring atoms depends on the 
oxygen coordination. From the symmetry considerations, 
an octahedral coordination provides better s-p overlap, 
Figure HJ compared to a tetrahedral one. Therefore, the 
energy dispersion of the conduction band is expected to 
be larger in the case of octahedral coordination. 

To elucidate the role of oxygen coordination in the 
charge transport properties of the oxides, we performed 
calculations for the oxygen deficient high-pressure a- 
G&2O3 phase with rhombohedral A^Oa-type structure 
(space group i?3c) and rock-salt ZnO phase. (In 2 03 with 
R3c (corundum) and 12 1 3 structures were not considered 
here because the cations (anions) have the same coordi- 
nation, i.e., octahedral (tetrahedral), as in bixbyite struc- 
ture.) 

In high-symmetry a-Ga2C>3 phase all Ga atoms have 
octahedral oxygen coordination - in contrast to fS-G^O^ 
with monoclinic C2/m structure which possesses both 
tetragonally and octahedrally coordinated Ga atoms. We 
find that in the oxygen deficient a-Ga2 03 with Vq con- 
centration of 1.3x 10 21 cm -3 , the nearest Ga and O atoms 
trap 53 % of the total charge density which is smaller 
compared to 61 % in (3-G&2O3 with similar defect con- 
centration of 0.8x 10 21 cm -3 . Consequently, the width 
of the single conduction band increases from 0.45 eV 
in /3-Ga 2 C>3 to 0.65 eV in a-Ga203 resulting in more 
than twice larger electron group velocity for the high- 
symmetry phase, Table HIl 



Similarly, we find that in oxygen deficient cubic ZnO 
with octahedral oxygen coordination of cations, the con- 
tributions from the atoms nearest to the defect decrease 
to 27 % of the total charge density and the contributions 
from atoms located at different distances from the oxy- 
gen defect differ by only 2 or 3 times. Accordingly, the 
conduction band width is larger (1.16 eV) that leads to a 
larger electron velocity near the Fermi level as compared 
to the ground-state hexagonal ZnO with the tetragonal 
oxygen coordination around Zn atoms (Table HI]) . 

Thus, we conclude that the high symmetry of octahe- 
dral oxygen coordination may improve conductivity ow- 
ing to a better s-p overlap which leads to a wider conduc- 
tion band and, hence, a larger electron velocity. A way 
to utilize the potential of an oxide is to ensure that the 
cations and anions possess a high local symmetry which 
can be attained via stabilization of a higher-symmetry 
phase or by choosing a multicomponent material with 
octahedral oxygen coordination for all constituents. 

III. 5 Role of cation's orbital contributions 

Despite the octahedral coordination of the cations in 
rhombohedral Ga203, the electron velocity in the oxy- 
gen deficient material still differ by almost a factor of 3 
from those in Iri2 03 and rock salt ZnO, Table HIl Again, 
the electron effective masses of stoichiometric oxides, cf. 
Table HI cannot explain the differences. We extend our 
investigation by comparing the cation's orbital contribu- 
tions to the conduction band. 

As mentioned in Section III.l, the cations and an- 
ions give nearly the same contributions at the bottom 
of the conduction band in all oxides investigated. Now 
we consider oxygen deficient materials and analyze the 
orbital contributions from cations, in particular, those 
that are nearest to the oxygen vacancy. (The contribu- 
tions from these atoms are considerable, judging by the 
charge build-up at the defect neighboring atoms which 
is observed for all oxides under investigation, see Section 
III.3 and Table HQ. Most importantly, we find that not 
only the s states but also the p and d states of the cations 
contribute to the bottom of the conduction band near the 
r point. This is expected from the kp theory, see for ex- 
ample [53J, which suggests that the states located both 
higher (conduction states) and lower (valence states) in 
energy with respect to a specified band (the conduction 
band minimum in our case) may contribute to that band. 
(Note that the sum in Eq. (1) runs over all bands except 
the conduction band). 

We compare the orbital contributions from the cations 
to the conduction band wave functions in the oxygen de- 
ficient In 2 03, rock salt ZnO and rhombohedral Ga 2 03. 
For the latter two oxides we use the high-symmetry 
phases where the cations have octahedral coordination. 
This will allow us to separate the role of local symmetry 
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FIG. 5. Band structure of undoped stoichiometric (a) In2C>3, (b) wurtzite ZnO, (c) /3-Ga203, (d) a-Ga203, and (e) InGaZn04 
as calculated within LDA. The empty p-states of cation(s) are marked with (+), (O) ° r (A) symbols to differentiate between 
the structurally and/or chemically distinct cations. Remarkably, the hybrid conduction band in InGaZn04 is formed well below 
the detrimental Ga p-states giving rise to a good charge transport. 



(the oxygen coordination) from the role of the electronic 
configuration of the cations in the observed differences 
in the transport properties. Values for the ground-state 
wurtzite ZnO and (3-G&2O3 are also given in Table [TTI for 
clarity. 

For In203 and rock salt ZnO we find that more than 
80% of the nearest cation contributions comes from their 
s-orbitals. Therefore, the conduction band originates pri- 
marily from the hybridization between the s-states of the 
cation and the antibonding p-states of the oxygen atoms. 
Clearly, the sp hybridization provides the most uniform 
charge distribution within the cell and, thus, facilitates 
good carrier transport. 

In marked contrast to In203 and rock salt ZnO, in 
rhombohcdral Ga203 the s-orbitals of the cations located 
near the oxygen defect contribute only 66%, while the 
contributions from the Ga p-states become considerable, 
21 %of the cations' total. This means that a substantial 
part of the charge introduced via the oxygen vacancy be- 
comes trapped on the p orbitals of the neighboring Ga 
atoms. The charge localization occurs due to the weak 
hybridization (small overlap) of these highly anisotropic 
p-orbitals with the p-orbitals of the oxygen atoms from 
the next coordination sphere. Low crystal symmetry re- 
duces this hybridization even further, as demonstrated 
by the results for /3-Ga203 and wurtzite ZnO, Table HH 

The origin of the suppressed s-orbital contributions in 
Ga203 and the resulting electron trapping near the defect 
is the proximity of the Ga p states to the conduction band 
bottom that gives rise to larger contributions from these 
p-orbitals. In Fig. \E\ we compare the calculated elec- 
tronic band structures for undoped stoichiometric oxides 
where the cation's empty p-states are highlighted. We 
find that the p-states in Ga203 are located more than 
1 eV closer to the bottom of the conduction band as 
compared to In203 and ZnO. Therefore, in contrast to 
indium and zinc oxides, the cation p-states in Ga203 may 



become available for the extra electrons - as we, indeed, 
found. 

Based on the above findings for single-cation TCOs, 
we can conclude that 

(i) the hybridization between the cation's s-orbitals 
and the p-orbitals of the neighboring oxygen atoms is 
crucial for achieving good charge transport; 

(ii) when the anisotropic p (or d) orbitals of the 
cation contribute considerably to the conduction band, 
the probability for the electron(s) to propagate further 
throughout the cell is reduced - owing to a weak hy- 
bridization between these p orbitals and the p orbitals 
of the neighboring oxygen atoms. As a result, the elec- 
tron velocity and, hence, conductivity are reduced. The 
degree of reduction depends on how close in energy the 
cation's p (or d) states are located with respect to the 
bottom of the conduction band. 

These findings explain the experimental observations 
that the conductivity of a multicomponent TCO de- 
creases with increasing Ga content. 



IV. Electronic properties of InGaZn04 

IV. 1 Undoped stoichiometric InGaZnOi 

Electronic band structure investigations of undoped 
stoichiometric InGaZn04 have been performed earlier for 
the crystalline 2j| 3(| and amorphous oxide 63]. Most 
importantly, it was found that despite the different band 
gap values of the constituent single-cation oxides, cf., Ta- 
ble |H all cations give comparable contributions to the 
conduction band bottom of the complex oxide. Because 
LDA may fail in description of the empty conduction 
states, the above finding of the hybrid nature of the con- 
duction band should be verified within sX-LDA. Here 
we performed such calculations for InGaZn04 and found 
that although the contributions (to the conduction band 
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at r point) from In states increase - from 23 % in LDA to 
33% in sX-LDA, those from Ga and Zn remain the same 
(15-16 % and 11 %, respectively) and so these states will 
be available for extra electrons once those are introduced 
into the host oxide. The change in the In contributions 
arises due to reduced contributions from 0(1) and 0(2) 
- namely, from 26 % in LDA to 23 % in sX-LDA and 
from 24 % in LDA to 17 % in sX-LDA, respectively. The 
reduction in the oxygen contributions is expected from 
a larger band gap in sX-LDA which we obtained to be 
equal to 3.29 eV, in good agreement with the experi- 
mental value of 3.5 eV (Til. Finally, the electron effec- 
tive mass remains isotropic in sX-LDA, Table HI confirm- 
ing that both structurally and chemically distinct layers, 
InOi.5 and GaZn02.5, will participate in charge trans- 
port upon a degenerate doping of this material. 




FIG. 6. Crystal structure of InGaZnCU. Only one of the 
three blocks which construct the conventional unit cell when 
stacked along the z direction, is shown. The cartesian z co- 
ordinates of the oxygen atoms are given on the left. 



IV. 2 Spatial distribution of oxygen defect in InGaZnOi 

First, we investigate the spatial distribution of the oxy- 
gen vacancy in the complex layered InGaZn04 by deter- 
mining the most energetically favorable location of the 
oxygen vacancy in the lattice. There are four struc- 
turally different sites for the defect, Fig. [6] and Table 
IIIII It can be located within the InOi.5 layer (this corre- 
sponds to a vacancy in the 0(1) structural type (38|) or 
GaZn02.5 layers (a vacancy in the 0(2) site). In addi- 
tion, because Ga and Zn atoms are distributed randomly 
in the GaZn02.5 layers [40], oxygen vacancy can have ei- 
ther Ga 3+ or Zn 2+ as an apical cation, i.e., the cation 
in the layer adjacent to the one where the defect is lo- 
cated. A comparison of the total energies of the oxygen 
deficient systems reveals that for all three charge states 
of Voj the defect prefers to be in the InOi.5 layer hav- 
ing Zn as its apical atom, Table IIIII We note that the 
energy difference between the oxygen vacancy located in 
InOi.5 and the one in GaZn02.5 layers increases signifi- 
cantly as the defect charge state decreases: it is only 56 
meV for uncompensated Vq, while for partially or com- 
pletely compensated defect the difference is an order of 
magnitude larger, Table IIIII 

An important observation is that the energetics of dif- 
ferent site locations of the oxygen defect in InGaZn04 
correlates with the formation energy of the constituting 
single-cation oxides, namely, -6.42 e\f /O2 for I112O3, - 
7.20 eV/0 2 for ZnO and -7.46 eV/0 2 for Ga 2 3 [55j]. 
The most preferred location of the vacancy (i.e., the 0(1) 
site with three In and one Zn neighbors) corresponds 
to the selection of the metal-oxygen bonds that would 
be easiest to break in order to create the defect. Al- 
though this conclusion must be verified for other systems, 
it may serve as initial guidance in defect studies of com- 
plex TCOs, both in crystalline and amorphous forms. 

From the formation energy calculations we find that 
the uncompensated oxygen vacancy, Vq , corresponds to 



TABLE III. Respective total energy, in meV, of InGaZn04 
system with uncompensated, partially compensated and fully 
compensated oxygen vacancies located at structurally differ- 
ent sites. 



Defect location V o V5 V 2 J 

Within the InOi.5 layer 

apical Ga +194 +138 +267 

apical Zn 

Within the GaZn02.5 layers 

apical Ga +110 +793 +1,319 

apical Zn +56 +536 +730 



the ground state of the defect - similar to the Ln^Oa, ZnO 
and Ga203 cases, Fig. [5J The LDA-derived energy re- 
quired to excite Vq to obtain Vq as well as the location 
of the oxygen defects with respect to the conduction band 
bottom in InGaZn04 are also similar to the single-cation 
oxides, Table |TTJ Hence, the vacancy can be photoex- 
cited to provide a metastable conductive state. Investi- 
gations of possible charge-compensating native defects 
(e.g., V/„ iG a,Zn or Zn Ga ) and defect complexes (e.g., 
(Zn<3 a Vo)) are required to further understand the lower 
observed conductivity in the layered InGaZn04 and in 
similar compounds from the series as compared to the 
single-cation TCOs. 

Finally, the obtained preferential distribution of the 
oxygen vacancies in InGaZn04 may explain the observed 
strong anisotropy of the electrical conductivity for the 
layered In203(ZnO) n compounds with an increased num- 
ber of ZnO layers, n >3 [3,[2l|,[2^]. A larger spatial sepa- 
ration between the InOi.5 layers due to extra ZnO layers 
may reduce the interaction between the oxygen defects 
which prefer to reside in the InOi.5 layers. This will lead 
to a smaller electron velocity along the z-direction. Wc 
note, however, that the anisotropy of the transport prop- 
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erties may not be well-defined in the as-grown samples 
since oxygen vacancy may occupy other, less favorable, 
locations. A higher but more anisotropic conductivity 
can be attained via a post-growth annealing with tem- 
peratures sufficient for the defect migration. Indeed, the 
observed significant decrease in the electrical conductiv- 
ity along the z axis with temperature which is reported 
to be almost twice smaller than the conductivity within 
the ab plane in the high temperature region [HI supports 
our conclusion. 



IV. 3 Electronic properties of oxygen deficient InGaZnOi 

Now we analyze the conduction band wave function in 
the oxygen deficient InGaZn04 - similar to the studies 
for single-cation oxides, see Section III. In accord with 
the hybrid nature of the bottom of the conduction band 
in the undoped stoichiometric InGaZnC>4, see [2{| and 
Section IV. 1, all atom types in the cell give compara- 
ble contributions to the conduction band wave function: 
In, Ga, Zn, 0(1) and 0(2) atoms contribute 25 %, 11 
%, 10 %, 32 % and 22 %, respectively. Moreover, for 
each type of atoms the contributions are fairly similar: 
the largest deviations are a factor of 2.5, 1.5, 2.0, 4.6 or 
1.6 for, respectively, In, Ga, Zn, 0(1) or 0(2) located 
at different distances from the oxygen defect. Therefore, 
the charge density in the conduction band of the oxygen 
deficient InGaZn04 is uniformly distributed throughout 
the cell - similar to 1^03 and ZnO cases and in contrast 
to /3-Ga203 where the defect's nearest and next near- 
est neighbor atoms give up to two orders of magnitude 
smaller contributions than the rest of the cations and an- 
ions in the cell. Accordingly, for multicomponent oxide 
we find that the energy dispersion of the conductive band 
and the electron velocities near the Fermi level are com- 
parable to those obtained for oxygen deficient L112O3 and 
superior to those for wurtzite ZnO, Table HH 

For InGaZn04 we find that the charge build-up at the 
defect's nearest cations and next nearest anions is 37 % 
which is notably lower than in Ga203 in either a or /3 
phase, cf., Table [TT1 It is important to point out that sim- 
ilar percentage is obtained for all locations of the oxygen 
defect, including the one where the vacancy has three Ga 
neighbors. 

Further, we analyze the relative orbital contributions 
from the defect neighboring cations to the conduction 
band wave function. We find that the s-orbitals of the In 
and Zn atoms nearest to the vacancy located in the InOi.5 
layer or of the Ga and Zn atoms nearest to the vacancy in 
the GaZn0 2 . 5 layer contribute 82% and 74% or 75% and 
71% of the cations total, respectively. The considerable 
increase in the s-contributions from the Ga atoms - as 
compared to the Ga203 phases, Tablellll- is attributed to 
the hybrid nature of the bottom of the conduction band 
in InGaZn0 4 . Due to the strong hybridization between 



the states of every cation with the states of its neighbor- 
ing oxygen atoms and because of the differences in the 
band gaps of the constituent oxides (3.1-3.6 eV for ZnO, 
3.5-3.7 eV for ln 2 3 and 4.4-4.9 eV for Ga 2 3 , found 
experimentally), the hybrid band is formed well below 
the Ga p-states, cf., Fig. [5fd). Qualitatively, this means 
that the interaction between the Ga s-states (conduction 
band) and the O p-states (valence band) become weaker 
in InGaZn04 (band gap 3.3 eV) as compared to Ga203 
(band gap ^4.6 eV), while the cation's p-states remain 
located at about the same energy, Fig. [5] As a result, the 
contributions to the conduction band wavefunction from 
these detrimental Ga p orbitals are notably smaller in 
InGaZn0 4 (14%) as compared to Ga 2 3 (33% or 21% in 
ft or a phase, respectively). Thus, the conduction band in 
complex InGaZn0 4 consists primarily from the cations s 
and the oxygen p states. This explains the obtained large 
electron velocity in this material, Table HU 

The above results highlight an important advantage of 
multicomponent oxides: mixing several oxides with no- 
tably different band gaps may help overcome the charge 
localization effects by tuning the resulting electronic 
properties via chemical composition of the complex ma- 
terial (and, to a lesser extent, via choosing/stabilizing a 
proper crystal structure). We must stress that both the 
relative content of the constituent oxides as well as their 
band gaps are critical parameters which control the re- 
sulting optical and transport properties. A large content 
of an oxide with a wider band gap compared to that of 
the other constituents may be appealing for the optical 
properties as well as for optimization of the material's 
work function desired for practical applications 64|. At 
the same time, however, as the band gap increases, the 
bottom of the hybrid conduction band will be effectively 
moved up - closer to the cation(s) p-states. This will in- 
crease the p-orbital contribution to the conduction band 
and limit the carrier mobility due to the increased prob- 
ability for the electrons to be trapped on the cation's 
p-states. 



V. Substitutional doping in InGaZnC>4 

Targeted doping with specific impurities which donate 
extra electrons upon aliovalent substitution represents 
an alternative route for the carrier generation in TCO 
hosts. Compared to the native point defects such as oxy- 
gen vacancy, substitutional doping is more appealing due 
to better stability of the doped samples to the environ- 
ment that is essential for practical applications. Among 
single-cation oxides, the most studied and commercially 
utilized are Sn-doped I112O3, Al-doped ZnO, Sb-doped 
Sn02 , F-doped ZnO and Sn02 . To the best of our knowl- 
edge, there were no reports on substitutional doping of 
InGaZn0 4 or other materials from the homologous series 
(in Ref. [ll, (ZnO) m In 2 3 was doped with Ca 2 +, Sr 2 + 
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TABLE IV. Properties of subtitutionally doped InGaZnCU. 
The dopant ionic radius, R, in A; the optimized sistances be- 
tween the dopant and its oxygen neighbors, D(M-O), in A; 
the formation energy of the defect, E/ orm , in eV; the electron 
velocities, v, calculated along the main crystallographic direc- 
tions in the vicinity of the Fermi level; the density of states 
at the Fermi level, N(Ef), in states/eV-cell; the fundamental 
band gap, E fundamental, in eV; the Burstein-Moss (BM) shift 
(the Fermi energy displacement measured from the bottom 
of the conduction band), in eV, for Ga, Sn, Ti and Zr doped 
InGaZnCU. The dopant concentration is 1.75xl0 21 cm -3 in 
all materials. For comparison, the ionic radii of In and Zn are 
0.94 A and 0.74 A, respectively. 



transition metal dopants provide an enhanced carrier mo- 
bility which leads to a pronounced increase in conductiv- 



ity with no changes in spectral transmittance [6JJ, |65j, |6 
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and Ba 2+ to improve the thermoelectric properties). 

One of the main challenges of substitutional doping in 
a multicomponent oxide is to find a proper dopant which 
will donate extra electron(s) and lead to degenerate dop- 
ing. The valence and the ionic size of both the dopant and 
all the cations in the cell should be taken into consider- 
ation. For example, doping of InGaZnC>4 with Al, which 
is known to be a good choice in the case of ZnO, is likely 
to result in Al 3+ (the ionic radius is 0.53 A) substituting 
Ga 3+ (0.61 A) instead of Zn 2 + (0.74 A) so that the ex- 
tra electrons are not generated and the material remains 
non-conducting. Moreover, non-uniform spatial distribu- 
tion of the dopants, their clustering or secondary phase 
formation as well as formation of defect complexes may 
become more pronounced in multicomponent oxides due 
to complexity and peculiarities in their crystal structure 
as compared to their single-cation counterparts. System- 
atic theoretical investigations may prove to be valuable 
to guide future experiments. 

Our goal is to investigate the structural, electronic 
and optical properties of doped InGaZnO.4 and also to 
compare them with those obtained for the oxygen defi- 
cient material. We begin by studying Sn 4+ and Ga 3+ 
substitutions. The former dopant is the most intuitive 
choice, while extra Ga may be considered as an intrinsic 
stoichiometry defect. In addition to these "traditional" 
dopants, we also considered Ti 4+ and Zr 4+ ions - fol- 
lowing recent experimental and theoretical findings that 



V. 1 Site location of Sn, Ga, Ti or Zr in InGaZnOi 

First, we determine which cation in the host struc- 
ture is most likely to be substituted by the particular 
dopant. For this, we compare the energies of the re- 
laxed structures where In or Zn is substituted by Sn 
or Ga; and In, Ga or Zn is substituted by Ti or Zr. 
Other configurations, e.g., Sn substituted on Ga sites, 
were not considered due to the large ionic radius mis- 
match. From the calculations of the defect formation 
energies, we find that Sn, Ti and Zr have a strong pref- 
erence to substitute In atoms whereas the substitution 
of Ga or Zn atoms corresponds to a higher energy by at 
least 1.4 eV, 1.3 eV and 3.6 eV, respectively. (The forma- 
tion energies are E/ or7n (Snz I j)=2.03 eV, E/ orm (Ti Ga )=- 
1.44 eV, E /orm (Ti z „)=-1.52 eV, E /orm (Zr Ga )=-0.29 
eV, E 

formi^Zn) — —0.33 eV. These should to be com- 
pared to the formation energies given in Table II VI) . 
The obtained preferential substitution can be understood 
based on the local symmetry arguments: ions with larger 
positive charge (M 4+ in our case) prefer to have more 
oxygen neighbors and so the octahedral oxygen coordi- 
nation (In sites) is more favorable than the tetrahedral 
one (Zn or Ga sites). Here we note that the obtained for- 
mation energies for doped InGaZn04 are similar to those 
in doped In203. For the latter, we find the formation en- 
ergy of the Sn substitution into In site to be equal to 0.17 
eV (which is in agreement with Ref. 58[), whereas for 
the Ti substitution E/ orm =-2.92 eV. 

In contrast, Ga shows a preference to substitute Zn 
atoms with the energy for In substitution being only 0.05 
eV higher (E/ orm (Ga/„)=-0.27 eV which should be com- 
pared to the value for Ga^„ , cf. , Table IIVI) . Because of 
the relatively small energy difference, one can expect that 
some of the extra Ga ions may substitute In atoms upon 
growth at elevated temperatures. This is a disadvantage 
since Ga has the same valence as In and so it will not 
contribute extra electrons upon the substitution. 

In Table IIVI the optimized distances between the 
dopant and its nearest oxygen neighbors are compared 
for doped InGaZn04 - along with the calculated elec- 
tron velocities and the optical properties to be discussed 
in Section V.3. 



InGaZnOi 



V.2 Electronic properties of 



We find that the electron velocity, Table IIVI is nearly 
isotropic in all doped materials - despite the non-uniform 
distribution of the dopants which prefer to reside within 
the InOi.5 layers. This is in accord with the hybrid 
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nature of the conduction band in InGaZnC>4 discussed 
above (see Section IV) and in our earlier work (29|. In- 
deed, in the Sn, Zr and Ga doped cases, the net con- 
tributions to the conduction band wave function from 
the atoms that belong to either In-0 or Ga-Zn-0 lay- 
ers are similar, namely, 58% or 42% in Sn case, 51% or 
49% in Ga and Zr cases, respectively. Furthermore, com- 
pared to the oxygen deficient material, see Section IV. 3, 
the charge density in the conduction band of the Sn, Zr 
and Ga doped InGaZn04 is more uniformly distributed 
throughout the cell since the largest deviations in the 
atomic contributions are only a factor of 1.2-1.6, 1.4-1.5, 
1.1-1.4, 1.2-1.7 and 1.1-1.6 for, respectively, In, Ga, Zn, 
0(1) and 0(2) located at different distances from the 
dopant. (Note that similar to the oxygen reduced com- 
plex oxides of the InGa03(ZnO)„ family, the isotropic be- 
havior may not be maintained in substitutionally doped 
materials with the increased number of ZnO layers (n) 
and anisotropic transport properties may appear after 
annealing.) 

The more uniform charge density distribution in the 
conduction band of Sn, Ga and Zr doped InGaZn04 leads 
to almost twice larger electron velocities (~8xl0 5 ra/s) 
than in the oxygen deficient material (4xl0 5 m/s) with 
the concentration of the partially compensated oxygen 
defect of 1.75xl0 21 cm -3 . (In this case, we compare 
the same defect concentration, with Sn or Vq giving one 
extra electron per impurity or per defect.) Similar re- 
sults are obtained for Sn-doped I112O3 where the elec- 
tron velocity is 1.4 times larger (9.17xl0 5 m/s for the 
dopant concentration of 1.93xl0 21 cm -3 ) as compared 
to the oxygen deficient indium oxide (6.42xl0 5 m/s for 
1.96xl0 21 cm -3 ). One should take into account the car- 
rier concentration in each case which is determined by 
the corresponding density of states (DOS) in the vicin- 
ity of the Fermi level. The DOS(Ei?) increases when 
the energy dispersion of the conduction band decreases 
which represents a more localized state. The conductiv- 
ity a depends on both, the electron group velocity and 
the density of states near Fermi: 



2c 2 

n 



^ \vk\\ 2 T k \S(E kx - E F ). 



(3) 



feA 



Here e is the electron charge, - the volume of the Bril- 
louin zone, k - the wave vector, A - the band index, v 
- the electron group velocity and Ef is the Fermi en- 
ergy. Based on our calculations, we can estimate the 
band-structure conductivity factor - the square of the 
electron velocity by the density of states, Eq. [3] The 
estimates suggest that the difference between oxygen de- 
ficient and Sn-doped oxides, InGaZn04 and In 2 03, are 
almost negligible. Therefore, the relaxation time for the 
Fermi surface electrons, t, will play crucial role in de- 
termining the resulting conductivity. While we cannot 
estimate t directly from the band structure calculations, 



below we speculate about the scattering mechanism asso- 
ciated with the Coulomb interaction between the ionized 
donor impurities and free electrons. (We do not consider 
the effects of grain boundaries, structural disorder and 
neutral impurities - which may have a significant effect 
on the transport prope rties.) 

According to 67H69J, the relaxation time depends on 
the concentration of impurity centers, Ni, which have 
charge Ze as well as the free electron density n. Because 
the charge state of the Sn 4+ substituting In 3+ ions, i.e., 
of the Sn defect, is the same as for the partially com- 
pensated oxygen vacancy, Vq, there is no apparent dif- 
ference in the relaxation times for the two carrier gen- 
eration mechanisms (Z—l and Ni=n for both). How- 
ever, we find that the structural relaxation is significantly 
stronger around the oxygen defect than near the substitu- 
tional dopant. For 1^03, the displacement of the nearest 
neighbor atoms from their original positions is 0.084 A 
and uniform around Sn ions, while it is much stronger and 
anisotropic around the partially compensated oxygen va- 
cancy with the displacements of 0.081 A, 0.112 A, 0.127 
A and 0.160 A for the four nearest neighbor In atoms. 
For InGaZn04, the nearest neighbor displacements are 
0.069-0.084 A in Sn-doped material and 0.097-0.331 A 
in the oxygen deficient oxide. Qualitatively, the larger 
atomic displacements around the defect are expected to 
give rise to significantly stronger scattering in the oxy- 
gen reduced oxides as compared to the substitutionally 
doped materials. Most importantly, a strong Coulomb 
attraction between the Vq defects and the free carriers 
and, hence, a short electron relaxation time, should be 
expected in all oxides considered in this work, since the 
partially compensated oxygen vacancy corresponds to a 
metastable state. 

Thus, proper doping with aliovalent substitutional im- 
purities is expected to lead to a notably higher electron 
mobility and conductivity than those observed in the oxy- 
gen reduced oxides. 



V.3 Role of the dopant's electronic configuration 

To elucidate the differences in the calculated electron 
velocities in Ga, Sn, Ti and Zr-doped InGaZn04, cf., Ta- 
ble IIV1 we compared the relative atomic contributions 
to the conduction band wave functions in these mate- 
rials. First we find that in the Ga, Sn and Zr cases, 
the anions give larger net contributions (64-66 %) com- 
pared to the cations (34-36 %), while in the case of the 
Ti doping the cations contribute twice as much (62 %) 
to the total charge in the cell. Moreover, in Ti-doped 
InGaZn04 about 64 % of the net cations ' contribution 
is due to the Ti d-states. In contrast, the states of Sn, 
Ga (primarily s-states) or Zr (primarily d-states) give 
only 18%, 5% or 7% of the total contributions from all 
cations in the cell. Therefore, the extra charge associated 
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with the aliovalent substitution is uniformly distributed 
throughout the cell in Sn, Ga and Zr-doped InGaZnGj, 
while it is mostly localized on the d-states of Ti [7(| re- 
sulting in the significantly lower electron velocity in the 
Ti-dopcd InGaZnO^. However, taking into account the 
increased density of states in the latter case, the band- 
structure conductivity factor, iV(e)u 2 (e), cf., Eq. [3l is 
approximately similar for all dopants investigated, with 
Ti showing the largest contribution to the band-structure 
conductivity, Table iTVl 

Further, we investigate how the dopant's electronic 
configuration and ionic radius affect the optical proper- 
ties of InGaZn04. In Table HVl the fundamental band 



gap and the Burstein-Moss (BM) shift [7l| are compared 
for doped materials. The fundamental band gap is mea- 
sured between the top of the valence band and the bot- 
tom of the conduction band (i.e., at T point), while the 
optical band gap in doped materials corresponds to the 
energy transitions from the top of the valence band to the 
Fermi level and, hence, is equal to the fundamental band 
gap plus the Fermi energy displacement (the BM shift). 
Note that, in contrast to oxygen deficient InGaZnC>4, the 
substitutionally doped material does not exhibit a second 
gap, as expected from hybridization between the dopant's 
states with the p states of neighboring oxygen atoms. 

We find that the BM shift is —0.4 eV smaller for the 
Ti-doped InGaZn04 as compared to Ga, Sn or Zr, Ta- 
ble IIV1 although Ti 4+ provides the same number of ex- 
tra electrons per substitution of In 3+ as Sn 4+ or Zr 4+ 
(or as the substitution of Zn 2+ with Ga 3+ ). This is ex- 
plained by the significant charge localization on the Ti 
d-states located in the vicinity of the Fermi level [70j |. 
The smaller BM shift should not affect the optical trans- 
parency within the visible range: because the band gap 
in InGaZn04 is large enough, 3.5 eV (Tflj], the energy of 
the transitions from the valence band is above the visible 
range (i.e., above 3.1 eV) even in the undoped oxide. 

Moreover, we argue that a smaller BM shift associated 
with the electron localization on the dopant d-states al- 
lows one to balance the optical and transport proper- 
ties more efficiently [6l|. In a conventional TCO, e.g., 
Sn-doped 1^03, an increased carrier concentration de- 
sired for an improved conductivity comes at a cost of 
an increased optical absorption because a pronounced 
BM shift reduces the energy of the transitions from the 
Fermi level, i.e., from the partially occupied conduc- 
tion band. In addition, plasma oscillations may affect 
the optical properties by reflecting the electromagnetic 
waves of frequency below that of the plasmon. In con- 
trast to the conventional (non-d element) doping, tran- 
sition metal dopants may allow one to introduce large 
carrier concentrations without compromising on the op- 
tical transparency. Indeed, smaller BM shift in Ti-doped 
InGaZn04 implies that the impurity concentration can 
be increased further (up to as high as 1.22-10 22 cm~ 3 , 
based on our electronic band structure estimates) before 



TABLE V. Optimized distances between the interstitial oxy- 
gen and its nearest cation neighbors, in A, and the respective 
total energy AE, in meV, for InGaZnCU systems with Oi lo- 
cated in various positions. The division into the groups is 
made based on the optimized cartesian z coordinate, see text 
and Figure [6] 
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reaching the same energy of the optical transitions from 
the Fermi level as with Ga, Sn or Zr doping at the dopant 
concentration of 1.75-10 21 cm -3 . We stress here that the 
above estimates are only to emphasize the trend between 
the s- and d-dopants. In a viable TCO, the impurity con- 
centration is limited by ionized impurity scattering on the 
electron donors as well as by the plasma frequency. 

Very recently, the experimental results for Mn doped 
InGaZn04 were reported [zH . It was found that the con- 
ductivity decreases while the carrier mobility increases 
with Mn content (up to 2.5 at.% of Mn). Such an in- 
triguing behavior may be associated with the onset of the 
long-range magnetic interactions between the Mn ions 
that leads to significantly longer relaxation time for the 
carriers of a particular spin [61]. At the same time, the 
density of states near the Fermi level, and thus, the car- 
rier concentration, decreases once the stable magnetic 
configuration is achieved. Accurate band structure cal- 
culations which go beyond LDA to correct the band gap 
value and the position of the localized Mn d-states with 
respect to band edges should be performed to verify this 
supposition. 



V-4 Charge compensation with interstitial oxygen 

It is known that high oxygen pressure can significantly 
reduce the conductivity of doped oxides due to the for- 
mation of neutral complexes between the ionized dopants 
and the interstitial oxygen atoms jsF 



73, [74 



Indeed, 

it was shown [6l| that interstitial oxygen in Mo-doped 
Ina03 may significantly affect the electronic, magnetic 
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and optical properties of the material. However, the 
formation energy of the Oi defect in undoped L12O3 is 
very high (we obtained Ey orm =8.88 eV for the oxygen- 
rich conditions that is in agreement with Ref. [58j). 
For a dopant-interstitial-oxygen defect complex (e.g., 
Mo"'Oj )* - given in Kroger- Vink notations) the energy 
of formation is lower by about 3 eV as compared to the 
interstitial oxygen alone [6lj . 

Given its compositional and structural complexity, 
InGaZnC>4 is a good test structure for investigation of 
the charge compensation mechanism via oxygen intersti- 
tial. Significantly, we find that the formation energy of 
the neutral O i defect in the multicomponent oxide is 4.77 
eV which is almost twice lower than that in ITI2O3. This 
suggests that the formation of dopant-O^ defect com- 
plexes in InGaZn04 is possible. While studies of particu- 
lar dopants and defect complexes in InGaZn04 is beyond 
the scope of this work, here we would like to understand 
the preferred distribution of the interstitial oxygen atoms 
in the complex oxide with chemically and structurally 
distinct layers. For this, we performed calculations for 
the undoped InGaZn04 with extra oxygen atom located 
in various structurally non-equivalent sites in the lattice. 
We find nine such positions for 0; which we divide into 
three groups according to their optimized Cartesian z 
coordinate: (i) the first three of them (cases 1, 2 and 3 
in Table fVf have z=3. 0-3.1 which is slightly above the 
framework oxygen atoms in the InOi.5 layer with z=1.9- 
2.2, cf., Fig. [6j (ii) another three with z=9.3 (cases 4 and 
5) and z=7.3 (case 6) are within the GaZn02.5 double 
layer where the framework oxygen atoms have z=6.2-6.5 
and z=9.9-10.2; (iii) the remaining three with z=4.7-5.5 
(cases 7, 8 and 9) reside in the "shared" space between 
the structurally and chemically different layers, and are 
about 2.5 times closer to the GaZn02.5 layer. 

Table |V] lists the respective total energies and the op- 
timized distances between the interstitial oxygen and its 
nearest cations for the structures with different location 
of the defect. As one can see, the interstitial oxygen 
prefers to reside within the InOi.5 layer having the dis- 
tance of ^2.1 A from one of the In atoms (cases 1, 2 
and 3). The total energy deviates insignificantly when 
the next nearest neighbor is either Ga or Zn. When the 
distance between the Oi and In atom increases (to 2.3, 
2.9 and 3.2 A) so does the total energy (by 127, 447 and 
1,490 meV, respectively). The energy decreases, how- 
ever, when the interstitial oxygen atoms reside within 
the GaZn02.5 layers. Nonetheless, the energy difference 
between the structures with located within the InOi.5 
or GaZn02.5 layers is at least 50 meV. Because the sub- 
stitutional dopants also have a preference to reside within 
the InOi.5 layers, as shown above, the formation of the 
charge compensated complexes is likely. Therefore, the 
ambient oxygen pressure must be carefully controlled to 
avoid the charge compensation and to ensure that the 
largest carrier concentration per substitutional dopant is 



attained. Using dopants which donate more than one 
electron per substitution, such as Mo 6+ [61J, represents 



an alternative route to mitigate the charge compensation 
and, thus, to achieve good electrical conductivity. 



VI. Conclusions 

The results of our comparative electronic band struc- 
ture investigations of In2 03, ZnO, Ga203 and InGaZn04 
allow us to draw several important conclusions: 

1) The electronic and optical properties of the oxides 
originate from the strong interaction between the metal 
s-orbitals and the p-orbitals of the neighboring oxygen 
atoms. Due to the Ms-Op overlap, the undoped stoichio- 
metric oxides exhibit large band gaps (3.4-4.9 eV) and 
small electron effective masses (0.28-0.35 m e ). The latter, 
however, cannot explain the differences in the observed 
electrical conductivities in degenerately doped oxides. 

2) Due to the three-dimensional Ms-Op network, all 
undoped stoichiometric oxides, including those with well- 
defined crystal lattice anisotropy, are capable to give rise 
to (nearly) isotropic conductivity. In layered InGaZn04, 
the electron effective masses are the same within and 
across the structurally and chemically distinct layers, as 
confirmed by our sX-LDA calculations. 

However, long-range structural anisotropy favors non- 
uniform (preferential) distribution of carrier donors. In 
/3-G8L2O3, oxygen vacancies are most likely to appear at 
the 0(1) sites which form zig-zag chains along the y axis 
leading to a twice larger electron velocity in this direc- 
tion. In InGaZn04, an oxygen vacancy, interstitial oxy- 
gen atom and substitutional Sn, Zr and Ti prefer to reside 
within the InOi.5 layer that may explain the observed 
temperature-dependent anisotropy of the electrical con- 
ductivity in InGa03(ZnO)„ and In203(ZnO)„ with an 
increased number of ZnO layers, n. 

3) Octahedral oxygen coordination provides the largest 
sp overlap, giving rise to a 1.6 and 2.5 times larger elec- 
tron velocity in hypothetical oxygen deficient rock salt 
ZnO and a-Ga203 as compared to wurtzite ZnO and 
/3-Ga203, respectively. Therefore, a way to improve con- 
ductivity and to utilize the potential of an oxide is to 
ensure that the cations and anions possess a high local 
symmetry which can be attained via stabilization of a 
higher-symmetry phase or by choosing a multicomponent 
material with, preferably, octahedral oxygen coordina- 
tion for all constituents. 

4) Most significantly, we find that the electronic config- 
uration of cation(s), in particular, the energy proximity 
of its empty p-states with respect to the conduction band 
minimum (CBM), plays the key role in determining the 
transport properties of oxygen deficient materials. In (3- 
Ga203, the Ga p states are energetically compatible with 
the s states of the host cations and, thus, available for the 
vacancy-induced electrons. As a result, in marked con- 
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trast to I112O3 and ZnO, the carriers become trapped on 
the p-states of the Ga atoms nearest to the defect - owing 
to the weak hybridization of these highly anisotropic Ga 
p-orbitals with the p-orbitals of their neighboring oxygen 
atoms. 

This finding explains why f3-G&2C>3 is not a viable TCO 
and also suggests that a large content of Ga203 in a mul- 
ticomponent TCO may lead to a suppressed conductivity 
- unless an alternative carrier generation mechanism is 
found. 

5) Due to the hybrid nature of the conduction band, 
multicomponent oxides offer a possibility to overcome the 
detrimental electron localization effects by tuning the re- 
sulting electronic - as well as optical - properties via 
chemical composition. 

6) In all oxides investigated in this work, neutral oxy- 
gen vacancy, V , results in a fully occupied, thus, non- 
conducting state located well below the oxide conduction 
band minimum. Because Vq also corresponds to the 
ground state of the defect, a conducting behavior may 
appear only if the vacancy is excited to Vq. In 1TI2O3, 
ZnO and InGaZnO.4, our LDA estimates suggest that the 
energy of excitation is within the visible range, i.e., the 
vacancy can be photo-excited, leading to a metastable 
conductive state. In P-G^O^, the excitation is highly 
unlikely because Vq and Vq prefer to be located at the 
different oxygen sites - of type 0(3) and 0(1), respec- 
tively. 

7) Aliovalent substitutional doping with next-column 
elements gives rise to twice larger electron velocity as 
compared to oxygen deficient materials. However, the 
carrier concentration is expected to be lower in the for- 
mer due to a smaller (by order(s) of magnitude) density 
of states in the vicinity of the Fermi level. This is asso- 
ciated with a nearly rigid-band shift of the Fermi level 
upon substitutional doping, whereas notable electron lo- 
calization near the oxygen vacancy leads to a higher DOS 
near the Fermi level in oxygen deficient materials. 

Our results suggest that the differences in the conduc- 
tivities in oxygen deficient and doped oxides stem from 
different carrier relaxation time. A shorter r in oxy- 
gen deficient oxides is expected from a significant atomic 
relaxation around the oxygen defect and a stronger 
Coulomb attraction between the oxygen vacancy Vq and 
the free carriers. 

8) Doping with specific transition metal elements rep- 
resents an efficient way to achieve high electrical conduc- 
tivity without substantial changes in the spectral trans- 
mittance. In addition, it provides a possibility to tune the 
electronic and optical properties of oxide via the forma- 
tion of charged defect complexes with interstitial oxygen 
atom. 

Thus, in-depth understanding of the microscopic prop- 
erties of the post-transition metal oxides provides signif- 
icant insights into the underlying phenomena in the con- 
ventional TCOs. Moreover, our systematic studies reveal 



several ways to knowledgeably manipulate the resulting 
properties - by tuning the chemical composition of com- 
plex hosts and carefully choosing proper carrier gener- 
ation mechanisms - to give rise to broader application 
range of these unique materials. 
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